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Abstract 

We apply periodic-orbit theory to calculate the integrated density of states N(k) from the peri- 
odic orbits of pseudointegrable polygon and barrier billiards. We show that the results agree so well 
with the results obtained from direct diagonalization of the Schrodinger equation, that about the 
first 100 eigenvalues can be obtained directly from the periodic-orbit calculations in good accuracy. 

PACS numbers: PACS numbers: 05.45.Mt, 



1 



The motion of a classical particle in a billiard system can show regular, chaotic or inter- 
mediate behavior, depending on the billiard geometry. A potential well of the same geometry 
as the corresponding classical billiard - a quantum billiard - reflects this behavior in the 
properties of its eigenvalues and eigenf unctions. A hallmark in the theory of chaotic systems 
is Gutzwiller's trace formula [l,Q|- It expresses the density of quantum mechanical eigen- 
states g(k) semiclassically by a weighted sum over all classical periodic orbits % of lengths 
li and thus represents an intrinsic link between the classical and the quantum mechani- 
cal properties of a given system. Since the implementation of Gutzwiller's trace formula, 
periodic-orbit theory has been a subject of permanent interest (for a recent review see j^l). 
However, in chaotic systems, the formula cannot be easily applied since the number N(£) 
of periodic orbits with lengths smaller than a given value £ increases exponentially with £, 
leading to a divergence of the trace formula. Therefore, cut-offs have to be compensated 
by sophisticated techniques so that the practical application of the trace formula up to now 
has been very limited and in most cases only reproduces the smoothed density of states and 
about the 20 — 40 lowest individual eigenvalues □ □ □ Only for the special case of the 
hyperbola billiard, where after a suited rearrangment of the orbits, most contributions were 
made to cancel, about 150 eigenvalues have been reproduced hi. 

In this paper, we concentrate on pseudointegrable billiards p[ , whose spectral proper- 



ties as e. g. t he 
billiards |loJ, DJJ 



evel statistics have been found intermediate between chaotic and integrable 
I12I [3]. Systems with rough boundaries may have chaotic or pseudointe- 



grable classical dynamics, so that the theoretical understanding of pseudointegrable systems 
is as important as the one of chaotic systems. One prominent application of quantum bil- 
liards is e.g. the band gap of Andreev billiards where the density of states is of major 
importance and where pseudointegrable billiards behave similar to chaotic billiards [l l 
Formulas equivalent to Gutzwiller's trace formula have also been established for regular [lj 
and for pseudointegrable billiards M . In these billiards, the number of periodic orbits smaller 

nnn 

than £ only increases as N(£) ~ £ 2 [ly, |12|, which diminishes the divergence problems. 
In this paper, we will show how the divergence problems can be overcome in these billiards 
and how the density of states and about the 100 lower eigenvalues can be calculated by 
periodic-orbit theory. 

Figure ^ shows some pseudointegrable geometries considered in this work together with 
some periodic orbits. Whereas rectangular systems are integrable, i.e. the motion of a 
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particle in a rectangular billiard shows regular classical dynamics and the equations of motion 
can be integrated, pseudointegrable billiards are polygons with a certain number of rational 

n 

angles <pi = riin/m i} with nj,mj e N and at least one n, > 1. Also barrier billiards [19| 
belong to this class (see Fig. ^d)), as a barrier can be considered as part of the boundary 
with an inner angle of 2ir. Pseudointegrable billiards are not integrable due to singularites 
arising at the salient corners and are classified by their genus number 

M J n- — 1 

G = l + — V^^. (1) 

Here, J is the number of angles and M is the least common multiple of the m^. In the 
geometries considered here, it is easy to see that every angle of value 37r/2 or 2n increases 
G by a value of 1, whereas the angles of ir/2 do not contribute. 

We are interested in the density of states g(k). It is well-known that in differential form 
g(k) reads 

g(k) = g (k)+g osc (k), (2) 

where k 2 = 2mE/h 2 and m is the mass of the quantum mechanical particle in the potential 
well. go(k) is a smooth term that can be obtained via the well-known Weyl formula j^j] 
from the geometrical properties of the system. It does not require the knowledge of the 
individual orbits. Hence, the calculation of g(k) reduces to the oscillating part g OS c(k) that 
is (for the billiards considered here) connected to the lengths and the areas of the periodic 
orbits. In both, integrable and pseudointegrable systems, the periodic orbits form families 
of fixed lengths £i, which means that the starting point of an orbit can be shifted to at least 
one direction along the boundary without changing its length (see Fig. Q). Accordingly, the 
trajectories of all orbit families cover a finite area Ai in phase space, where the index i counts 
the different families. Gutzwillers trace formula has also been extended to pseudointegrable 
billiards J^l Q]. For the billiards considered here, where all orbits have an even number 
of reflections at the boundary walls, the boundary conditions do not play a role and the 
formula reads 

9osc{k) = J^Y,jk C °^ Ml ~^ (3) 
and applies also for integrable billiards. Second-order contributions coming e.g. from diffrac- 
tive orbits have been neglected and the sum is carried out over all primitive (non-repeated) 
orbit families i and over its repetitions with multiple lenghts. Since the boundary conditions 
do not enter into Eq. ©, we concentrate in the following on Dirichlet boundary conditions. 



Even though the number of orbits N(£) below a given length £ increases only quadratically 
with £, one can easily show that also the trace formula for pseudointegrable systems, Eq. @, 
diverges and therefore could not be used so far to calculate the density of states. However, 
since the divergence is weak, we can use a simple trick to achieve convergence, namely by 
considering the fluctuations of the integrated density of states, 

N osc (k) = / g osc (k) dk = 



I \ I k 

" ' 7^3 [sin(^fe) - cos(£ik)} 




where FrS(x) and FrC(x) are the Fresnel sine and cosine integrals, respectively, which can be 
evaluated numerically. Replacing the sum over % by an integral over g(£) d£ and introducing 
the orbit density g(£) = dN(£)/d£ ~ £, one can verify easily that the factors of £^ 3 ^ 2 and 
£^ 2 ensure the convergence of Eq. (jlj). 

We first test the formula on the rectangular billiard where both, the orbits and the 
quantum mechanical eigenvalues k 2 are known exactly, k 2 = n 2 \v 2 j 'L 2 + v 2 / L 2 ) with 
positive integers z/j and Uj, and the side lengths L x and L y of the rectangle. The orbit 
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lengths are £ Vi ,v j = 2 [(viL x ) 2 + (ujL y ) 2 ] and the areas are Ai = 2A for the neutral orbit 
families (the simplest orbit families that bounce between two parallel walls 3|) and 4A for 
all other families, where A is the geometrical area of the system. In Figs. Efa,b) we compare 
the integrated density of states N PO {k) = N (k) + N osc (k) (straight lines) calculated from 
Eq. with the corresponding density of states N EY (k) (circles), which has been obtained 
from the exact eigenvalues. No(k) is taken from Weyl's formula (for Dirichlet boundary 
conditions) , 

where A is the area, T the boundary length of the billiard and the sum runs over all corners 
of angles One can see in Fig.^lthat N PO (k) gives the expected staircase function and the 
agreement to N EY (k) is excellent. As shown in the upper curve of Fig. Efb), it even allows 
to obtain the eigenvalues directly from the steps in N p (k). Each eigenvalue is positioned 
at one of the steps of N PO (k), which we verified for the first 1500 eigenvalues (until k 2 = 1). 

In pseudointegrable billiards, the areas are different for the different families and normally 
much smaller than in integrable billiards, while the number of periodic orbits is larger. 



In Fig. 1(b), it is demonstrated how the periodic orbits can be labeled according to their 
numbers of transversals of the different segments: the system has two x- and two y-segments 
and the orbit shown here can be labelled as (2, 6, 2, 0), where the numbers design the number 
of transversals of the segments x±, X2, y\ and y<i (see also Ref. ^|). The orbit length 
and angle can be calculated from this information, whereas the area has to be calculated 
numerically. However, unfortunately not all combinations of integer transversals exist in 
the pseudointegrable systems, because many hypothetical orbits are pruned by the shielding 
of the corners. So, it can be seen easily in Fig. 1(b) that an orbit (2,0,2,0) would not 
be possible in this geometry (however, this orbit exists in the geometry of Fig. 1(a)). Even 
though the numerical calculations can be done in high precision for up to about 50000 orbits 



in reasonable computation time [22j, we first want to investigate the stability of the results 
against numerical errors. To this end, we compared for the pseudointegrable billiards the 
orbits found in forward and backward direction. The maximum errors of the areas of the 
first 40000 periodic orbits (about 300 reflexions at the system walls including repetitions) 
are around 0.1 % and we found no hints that orbits might be lost. Accordingly, we have first 
tested the robustness of the results by disturbing the orbits of the rectangular system by 
errors taken from a narrow Gaussian distribution of a = 0.05 and a maximum error of 0.1% 
(see Fig. Efb) shifted down by a value of 5 in the lower curve (dotted line)). One can see 
that on the average, N PO (k) still agrees very well with N EV (k), but that the shape of the 
staircase is smeared out already by these quite small errors, so that the eigenvalues can no 
longer be determined from the steps in N PO (k). The loss of about 1% of the orbits (as well 
as their repetitions) would be less disturbing and would only lead to very slight deviations 
in the step function. We also checked if the number of calculated orbits is large enough and 
found that for all considered systems, the results are stable beyond the first 10 4 orbits. 

Finally, we turn to the pseudointegrable billiards. As before, we apply Eq. (JU) to obtain 
the density of states N (k) from periodic-orbit calculations. For comparison with the 
quantum mechanical eigenvalues, we solve the Schrodinger equation for the potential wells 
of Fig. [U in its discretized version, 

£ (^j") -*»(*, J)) = "(^) *»(»*. J). (6) 

V ,j' (neighbors) \ / 

where the indices refer to points of a square lattice with v lattice points per unit length. 
The sum over (i r ,f) runs over all nearest neighbors of Equation © depends on the 



usual Schrodinger equation, A^f n (i,j) = —k 2 ^/ n (i, j), by a second-order Taylor expansion 
of the left-hand side up to the quadratic term. The errors due to the discretization are the 
higher orders of v and thus decrease with increasing u, while the eigenvalues are transformed 
via k 2 —>■ ik n jv) 2 - We used a resolution of v — 4, where the errors arising from the 
discretization (as calculated for the rectangle) are smaller than 0.03 percent for the first 100 
eigenvalues. Equation © describes a matrix problem that we diagonalized by the Lanczos 
algorithm, yielding the numerical Lanczos eigenvalues k\ and the corresponding density of 
states N ' (k). In Fig. 01 we compare for the polygone billiards of Figs. H|a-c) N PO (k) 
with N EV,L (k). First, in Fig. 0(a), we show a larger part of the energy spectrum, where 
we can see that both, N PO (k) (straight lines) and N EV,L (k) (symbols) agree again very 
well. Even though the steps in N PO (k) are smeared out by the numerical inaccurencies as 
described above, we can use them to obtain the individual eigenvalues from N p (k). To 
this end, we fitted N p (k) by a least square fit to a step function with constant integer 
values of the step heights. The positions kf of the steps were chosen by minimizing the 
quadratic deviation to N PO (k). The step functions obtained this way, N p t °(k), are shown 
in Fig. 0(b) (solid lines). We can see that indeed, for about the first 100 eigenvalues, the 
agreement of the Lanczos eigenvalues with this fitted step function N p t °(k) is quite good. 
Nearly all eigenvalues are located right at the steps of N p t °(k), and we can obtain at least 
the first 100 eigenvalues k 2 PO by periodic-orbit theory. In Tab. HJ we show as an example the 
91st to the 100th eigenvalue. We calculated the number N m of mismatches among the first 
100 levels, i.e. the number of cases where the ith periodic-orbit eigenvalue lies closer to the 
(i + l)th or (i — l)th Lanczos value than to the ith one. N m is also given in the table and lies 
around 20 percent. Naturally, the mismatches occur at energies, where the level distance is 
particularly small and N m is therefore highest for the system of Fig. ^a), where the density 
of states increases fastest and the level distances are thus smallest. 

As last example, we consider the barrier billiard of Fig. E^d) for three different heights 
of the barrier, h — 10, h — 50 and h = 100. In this case, the discretization of the lattice 
is a cruder approximation than before, since the barriers which should be of thickness zero, 
always occupy one grid point. In Fig.0J we can see that at least for the systems with barrier 
heights h = 10 and h = 50, the agreement between the periodic-orbit and the Lanczos 
results is again very good with mismatches even below 20 percent (see Tab. |HJ). Only for 
the billiard with the largest barrier height, the mismatches are larger which is, however, most 
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probably due to the discretization procedure and not to the periodic-orbit calculations. A 
comparison to the eigenvalues calculated by some other procedures will be interesting. 

Finally, we want to compare our results to some other methods that were used in the 
past to determine eigenvalues of systems with chaotic classical dynamics by periodic-orbit 
theory. The first methed used e.g. in \^ is a "Gaussian smearing" of g(k). This means that 
convergence of the trace formula can be achieved by multiplying each term of Eq. (J3j) with the 
additional factor of exp(— lfe 2 /2) (where e must be small). As a consequence, the delta peaks 
of g(k) are transformed into Gaussian functions of shapes ~ exp(— k 2 e~ 2 /2). We tested this 
method also on our systems and found that it works very well for the integrable rectangular 
systems, but not for the pseudointegrable systems. Similar as in our method, the small errors 
in the orbit areas lead to a large noise that is still increased by the "Gaussian smearing". 
Contrary to our method, where we were able to apply a clearly defined fit procedure to 
eliminate the noise from the steps in N(k), it is not possible to find the Gaussian functions 
in g(k) by a simple rule. Another method used in 3| calculates the eigenvalues from the 
zeros of the so-called dynamical zeta function that contains all orbit information. In very 
special cases, where every large primitive orbit can be decomposed into series of few small 
orbits, the dynamical zeta function can be calculated very easily. However, this condition 
is only fulfilled in rare cases and demands as minimal conditions that the orbits can be 
labelled into fundamental building blocks, where each combination of the blocks exists. We 
have seen that in the case of pseudointegrable billiards, the periodic orbits appear in a 
very unsystematic way and many hypothetical orbits are pruned, so that we think that this 
method is not helpful in our case. 

In summary, we have shown how the convergence problems of the trace formula can be 
overcome in systems, where the number of periodic orbits below a given length I increases 
at most quadratically with £, e.g. for integrable and for pseudointegrable billiards. We 
have shown that the integrated density of states can be reproduced in very good accuracy 
for several hundred eigenvalues. The calculations are very sensible to numerical errors, so 
that already error bars of about 0.1% destroy the shape of the step function of N PO (k) 
such that the steps are smeared out. In integrable billiards, where the orbits are known 
exactly, the eigenvalues can be found directly from the steps in N PO (k) (which we tested 
for the first 1500 eigenvalues). In pseudointegrable billiards, even if the curves follow all 
fluctuations of the spectra very well, a fit technique has to be used in order to find the 
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first 100 individual eigenvalues. The results are very promising and show that the quantum 
mechanical density of states can indeed be gained, using as input solely the classical periodic 
orbits of the pseudointegrable billiards. It will be very interesting to apply this procedure to 
real mesoscopic structures, where the change of the density of states as response to a change 
in geometry is of great importance. 

We would like to thank A. Bunde for a careful reading of the manuscript and valuable 
remarks. 
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FIG. 1: Shapes of the pseudointegrable billiards considered in this work. (a,b) L-shaped billiards 
(polygon billiards of genus number G = 2), (c) a polygon billiard of G = 3 and (d) the barrier 
billiard (G = 2), which is calculated for different heights h of the barrier. One periodic orbit for 
each geometry is also shown (dashed lines). In (b) the different segments are shown that can be 
transversed by the periodic orbits. 
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FIG. 2: (a,b) Density of states N PO (k) calculated from the periodic orbits (solid lines) and N EV (k) 
calculated from the exact eigenvalues (circles) of a rectangular system of side lengths L x = 101 
and L y = 198. In (b), the upper curve compares the steps in N PO (k) to the eigenvalues. In the 
lower curves (shifted down by a value of 5), the areas of the orbits are disturbed by errors of up to 
0.1% and it can be seen that the steps in N PO (k) (dotted line) are strongly disturbed already by 
these small error bars. 
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FIG. 3: (a,b) Density of states N PO (k) calculated from the periodic orbits (solid lines) and 
N EV,L (k) from the Lanczos eigenvalues for the L-shaped systems of Fig.^a) (circles) and Fig.^b) 
(squares) and the system of genus number G = 3 of Fig. ^c) (diamonds). In (b), the fitted step 
function of N^ t °(k) is compared to the Lanczos eigenvalues. For a better overview the data of the 
geometries of Figs. ^b,c) have been shifted upwards by values of 150 in (a) and by 25 in (b). 
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FIG. 4: (a,b) Density of states N PO (k) calculated from the periodic orbits (solid lines) and 
N EV,L (k) from the Lanczos eigenvalues for the barrier billiards with barrier heights h = 10 (di- 
amonds), h = 50 (squares) and h = 100 (circles). In (b), the fitted step function of N^ t °(k) is 
compared to the Lanczos eigenvalues. For a better overview the data in both, (a) and (b) have 
been shifted upwards by values of 5 (h = 50) and by 10 (h = 10). 
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TABLE I: Table of the 91st to the 100th eigenvalue kp calculated by the trace formula compared 
to kl calculated by the Lanczos algorithm for the polygon pseudointegrable systems. N m is the 
number of mismatches during the first 100 values, i.e. the number of cases where the ith periodic- 
orbit eigenvalue lies closer to the (i + l)th or (i — l)th Lanczos value than to the ith one. 
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TABLE II: Table of the 91st to the 100th eigenvalue k PO calculated by the trace formula compared 
to kl calculated by the Lanczos algorithm for the pseudointegrable barrier systems. N m is the 
number of mismatches during the first 100 values. 
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